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Abstract 



A systematic investigation on vapor-liquid equilibria (VLE) of dipolar and quadrupolar 
fluids is carried out by molecular simulation to develop a new Helmholtz energy contri- 
bution for equations of state (EOS). Twelve two-center Lennard- Jones plus point dipole 
and point quadrupole model fluids (2CLJDQ) are studied for different reduced dipolar 
moments //*^ = 6, 12, reduced quadrupolar moments Q*^ = 2, 4 and reduced elongations 
L* = 0, 0.505, 1. Temperatures cover a wide range from about 55 % to about 95 % of 
the critical temperature of each fluid. The NpT -\- test particle method is used for the 
calculation of vapor pressure, saturated densities and saturated enthalpies. Critical data 
and the acentric factor are obtained from flts to the simulation data. On the basis of this 
data, an EOS contribution for the dipole-quadrupole cross-interactions of non-spherical 
molecules is developed. The expression is based on a third-order perturbation theory 
and the model constants are adjusted to the present 2CLJDQ simulation results. When 
applied to mixtures, the model is found to be in excellent agreement to results from simu- 
lation and experiment. The new EOS contribution is also compatible with segment-based 
EOS, hke the various forms of the Statistical Associating Fluid Theory (SAFT) EOS. 
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1 Introduction 



Knowledge of thermodynamic properties, and in particular vapor-liquid equilibria (VLE), 
is important for many problems in science and engineering. Among the state-of-the-art 
thermodynamic models, molecular based approaches have the highest potential to yield 
significant improvements compared to existing phenomenological approaches, especially 
in terms of predictive power. 

For the two-center Lennard- Jones plus point dipole (2CLJD) fluid and its quadrupolar 
pendant (2CLJQ), systematic studies of VLE were carried out in previous work |T1 12]. 
The VLE results were correlated as functions of the model parameters, which were used 
to determine the model parameters for 78 real fluids [H H] . It has been shown that these 
molecular models can successfully be applied for the description of VLE for binary and 
ternary mixtures OElEj. The simulation data from the two systematic studies pQI2] were 
subsequently used to construct equation of state (EOS) contributions due to dipolar [8J 
and quadrupolar [9] interactions. 

As numerous real fluids, e.g., carbon monoxide or refrigerants like R115 (CF3-CF2CI) 
are both dipolar and quadrupolar, it is valuable to study multipolar model fluids. In 
case of mixtures, containing dipolar components and quadrupolar components, the polar 
cross-interaction also plays a significant role. However, only few studies on VLE of fluids 
that are both dipolar and quadrupolar are available in the literature. Without reference to 
real fluids, Dubey and O'Shea studied VLE of one-center Lennard- Jones plus dipole 
and quadrupole model fluids. 

Molecular theories for multipolar compounds were developed by Stell et al. [HI \n\ . 
based on a perturbation theory around a non-polar reference fluid. A hard-sphere refer- 
ence was considered by Rushbrooke et al. p3] and later by Henderson et al. [H]. The 
perturbation theory for mixtures was worked out by Gubbins and Twu [IS] consider- 
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ing point-electrostatic multipoles with the spherical Lennard- Jones reference fluid. The 
perturbation terms were initially parameterized to molecular dynamics results and this 
parameterization was subsequently improved by Luckas et al. The perturbation 

terms were in detail reported by Moser et al. [17] and Shukla et al. ^18j and good results 
were found upon applying that theory to real mixtures [19l EH] . Boublik studied the ef- 
fect of the non-spherical molecular shape on multipolar interactions, considering the pair 
correlation function of a Gaussian overlap fluid [21] . 

An equation of state for the 2CLJD and the 2CLJQ fluid were developed by Saager 
and Fischer based on molecular simulations. A fixed elongation {L* = 0.505) was thereby 
considered and a dipole-dipole term and a quadrupole-quadrupole term were obtained 
by fitting empirical expressions to the simulation data [221 EH] . Dipole-quadrupole cross- 
interactions were subsequently treated with an effective one-fluid dipole moment which 
contains a contribution due to the quadrupole moments of the mixture [21] . Similarily an 
effective one-fluid quadrupole moment was defined where the dipole moments of compo- 
nents in the mixture contributed according to an empirical relation that was parameterized 
to simultion data. 

Due to the success of the Statistical Associating Fluid Theory (SAFT) EOS [25112611271 
[2H1 in various variants [301 [311 E21 E31 EH ES] , the multipolar terms developed earlier 
for spherical fluids, were recently applied in combination with different SAFT models 
[36l[371[38l[39lll0llillll2lll3llll]. Zhao and McCabe used an integral-equation approach 
with the mean spherical approximation (MSA) closure in combination with the SAFT-VR 
EOS [45j. As structural information is available from the MSA term, the orientation of 
dipoles can be accounted for. 

The authors of the present study have recently developed expressions that account 
for the non-spherical shape of fluids and found good results also for mixtures of 2CLJQ 
and 2CLJD fluids as well as for mixtures of real substances [8l|9]. Furthermore, Kleiner 
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and Gross [IB] accounted for the polarizability of fluids and for the appropriate induction 
effects of dipoles in applying a renormalization scheme proposed by Wertheim [171 H8l HHj 
l3Uj . A combination of the perturbation theory of Moser et al. [T7] and of Shukla et al. 
[18] . where the full vectorial and tensorial information of multipoles is utilized, with the 
approach of our earlier work [HI [9] is presented by Leonhard et al. [5l] . 

In this work, VLE of two-center Lennard- Jones model fluids which are both dipolar and 
quadrupolar (2CLJDQ) are systematically studied over a wide range of model parameters 
and temperatures. In reduced units, the 2CLJDQ model class has three molecular pa- 
rameters that can be varied: elongation, dipole moment and quadrupole moment. Twelve 
different 2CLJDQ model fluids are covered here. Using these results, an EOS contribu- 
tion for the dipole-quadrupole cross-interactions of non-spherical molecules is developed 
based on a third-order perturbation theory. Model constants are thereby adjusted to the 
simulation results of the 2CLJDQ fluid. Finally, the EOS contribution is compared to 
experimental and simulation data for five real binary mixtures. 

2 Molecular model 

The two-center Lennard- Jones plus point dipole and point quadrupole fluids (2CLJDQ) 
is here investigated. The pair potential of this model class is composed of two identical 
Lennard- Jones sites a distance L apart (2CLJ) plus a point dipole of moment /i and a 
point quadrupole of moment Q placed in the geometric center of the molecule. Both 
polarities are aligned along the molecular axis. The full potential writes as 

M2CLjDQ(»'ii, ^^u ^j, L, fi, Q, a, e) = U2CLj{rij, uji, uj, L, a, e) + UBD{rij,u;i, ujj, /i) 
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where 

2 2 
a=l 6=1 

is the Lennard- Jones part. Herein, is the center-center distance vector of two molecules 
i and j, Tab is one of the four Lennard- Jones site-site distances, a counts the two sites of 
molecule i, b counts those of molecule j. The vectors uji and loj represent the orientations 
of the two molecules. There are three polar contributions, cf. Allen and Tildesley [52] : 
Firstly, the interaction between dipoles is present 

UBD{rij, uji, u}j,fi) = ^ {siSjC - 2ciCj) , (3) 

secondly, the interaction between quadrupoles 

3 2 
'^QQif ij, ^i, ^j, Q) = -7- ^ [1-5 [c^i + cj) - ISq^cJ + 2 {siSjC - Aacjf] , (4) 

and finally, the mutual interaction between dipoles and quadrupoles 

3 

u-DQirij,uJi,ujj,fi,Q) = {d - Cj) {1 + SciCj - 2siSjc) , (5) 

where = cosd^, = sinOk, and c = cos^jj. di is the angle between the axis of the 
molecule i and the center-center connection line and (pij is the azimuthal angle between 
the axis of molecule i and j. The Lennard- Jones parameters a and e represent size and 
energy, respectively. 

Among the different spatial arrangements of the four charges in a quadrupole, in 
2CLJDQ models they are arranged along the molecular axis in the symmetric sequence 
-|-, — , — , -|- or, having the same energetic effect in pure quadrupolar fluids, — , +, +, — . 
The point quadrupole interaction, cf. Eq. (jl]), is a large distance approximation, reducing 
the number of parameters related to the quadrupole to one, namely the quadrupolar 
moment Q. The dipole is approximated in the same sense. 
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In case of elongated fluids, for very small intermolecular distances \rij\, the positive 
Lennard- Jones part U2clj of the full potential cannot outweigh the divergence to — oo of 
the polar part mdd+wqq+mdq. This divergence of 'U2Cljdq leads to infinite Boltzmann 
factors, i.e. non-existence of the configurational integral. During molecular dynamics 
phase space sampling within the pressure range in question, this artifact of the 2CLJDQ 
potential causes no problem as intermolecular center-center distances are very improbable 
to fall below critical values. However, during Monte-Carlo simulation or the calculation 
of entropic properties by test particle insertion [53], critical intermolecular center-center 
distances might occur. To avoid this problem, following Moller and Fischer [SI], a hard 
sphere of diameter OAa was placed directly on the polar sites as a shield for critical 
configurations. 

The parameters a and e were used for the reduction of all thermodynamic properties 
as well as the model parameters: T* = Tks/e, p* = pcr^/e, p* = pa^, h* = h/e, L* = 
L/a, /i*^ = /i^/ (ecr^) and Q*^ = Q"^/ (ecr^). The model parameters were varied in this 
investigation: L* = 0, 0.505 and 1, /i*^ = 6 and 12 as well as Q*'^ = 2 and 4. Combining 
these values leads to a set of twelve model fluids that were investigated here. 

To achieve a monotonous transition from L* > to L* = 0, spherical fluids were 
treated as 2CLJ fluids with L* = 0. This leads to a site superposition that is not present 
in the one-center Lennard- Jones case. Therefore, in reduced units, temperature, pressure 
and enthalpy are fourfold and the squared polar moments twofold of the corresponding 
values when only one Lennard- Jones site is present. Densities are not affected. 

3 Molecular simulation method 

Throughout all pure fluid VLE simulations, the NpT + test particle method by Moller 
and Fischer [551 156] was used. The chemical potential was calculated in the NpT ensemble 
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by Widom's method [SB] . 

In all simulations a center-center cut-off radius = Scr was used. Outside the cut-off 
sphere, the fluid was assumed to have no preferential relative orientation of the molecules, 
i.e. for the calculation of the Lennard- Jones long range corrections, orientational averaging 
was done with equally weighted relative orientations as proposed by Lustig Long 
range corrections for the dipolar interaction were calculated with the reaction field method 
[HH1I22!, where the relative permittivity eg was set to infinity. The quadrupolar interaction 
needs no long range correction as it disappears by orientational averaging. The same holds 
for the mixed polar interaction between dipoles and quadrupoles, cf. Weingerl and Fischer 

m- 

Configuration space sampling was done with = 864 particles for both liquid and 
vapor simulations in the NpT ensemble with Andersons barostat [59]. The reduced in- 
tegration time step was set to l^t^fmftj o = 0.0015 and the reduced membrane mass 
parameter of the barostat was set to 2 ■ 10~^ for liquid and to 10^^ for vapor simulations. 

Starting from a face centered lattice arrangement, every run was equilibrated over 
10, 000 time steps. Data production was performed over 100, 000 time steps. For vapor 
simulations with L* = 1, especially those at low temperatures, the Monte-Carlo method 
was used to achieve equilibrium within an acceptable time range. The number of Monte- 
Carlo loops was generally chosen to be the same as the number of MD time steps. One 
Monte Carlo loop is defined here as N trial translations, (2/3) trial rotations, and 
one trial volume change. At each production time step 2A^ test particles were inserted 
in the liquid phase, and A^ test particles in the vapor phase to calculate the chemical 
potential. To get a better accuracy of the chemical potential, iN test particles were used 
in simulations in the liquid phase at low temperatures. 

In some cases, where a highly dense and strongly polar liquid phase was present. 
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the more elaborate gradual insertion scheme had to be employed to obtain the chemical 
potential sufficiently accurate. The gradual insertion method is an expanded ensemble 
method [20] based on the Monte Carlo technique. Here, the version proposed by Nezbeda 
and Kolafa [61] was used, in a form that was extended to the NpT ensemble [62]. In 
comparison to Widom's test particle method, where whole molecules are inserted in the 
fluid, gradual insertion introduces one fluctuating molecule, that undergoes changes in a 
predefined set of discrete states of coupling with the other molecules of the fluid. Pref- 
erential sampling is done in the vicinity of the fluctuating particle. This concept leads 
to considerably improved accuracy of the residual chemical potential. Gradual insertion 
simulations were performed with N = 864 particles in the liquid phase. Starting from 
a face-centered lattice arrangement every simulation run was given 5, 000 Monte Carlo 
loops to equilibrate. Data production was performed over 100, 000 Monte Carlo loops. 
Further simulation parameters for runs with gradual insertion were taken from Vrabec et 

ai. [sg. 

VLE data were determined for temperatures of about 55 % to 95 % of T*{fi*^, Q*'^, L*). 
In the whole temperature range all thermodynamic properties of both phases were ob- 
tained by simulation. 

4 Simulation results of VLE data 

Table [T] reports the VLE data of the regarded twelve model fluids. Vapor pressure p* , 
saturated liquid density p'*, saturated vapor density p"*, residual saturated liquid enthalpy 
^/res* g^j^^ residual saturated vapor enthalpy h"'^^^* are presented, where h^^^* = h*{T,p) — 
h^'^{T). Statistical uncertainties were determined with the method of Fincham et al. [63j 
and the error propagation law. 

Figs. [H to [H] illustrate the strong influence of both elongation and polar moments on 
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the VLE data. Elongation and polar moments strongly influence saturated density and 
the vapor pressure curves. As can be expected, all three critical properties in Figs. [1] to 
[HI i.e. temperature, density and pressure, decrease with increasing elongation for a given 
set of polar moments and thus the phase envelopes migrate accordingly. For a constant 
elongation, they usually increase with increasing polar moments. 

To determine the critical data quantitatively, the method of Lotfi et al. [M] was used. 
The saturated density-temperature dependence near the critical point is well described by 
p* ~ (T* — T*)^^^, as given by Guggenheim [651 [66] • Therefore, the following correlations 
were fitted to simulated vapor pressure and saturated densities 

Inpl = ci + || + ^, (6) 
p'* = pI + Ci ■ (T; - T*y/^ + C'^ ■ {T* - T*) + C'^ ■ {T* - T*)3/2, (7) 

p"* = pi-Ci- {T* - T*)i/=^ + c'J, ■ (t; - T*) + c'^ ■ (t; - t*)^/^. (s) 

The simultaneous fit of saturated densities, cf. Eqs. ([7j) and ((HI), yields p* and T*. Inserting 
T* into the vapor pressure fit, cf. Eq. ([H]), gives the critical pressure p*. These fits can be 
seen in Figs. [UtolHl The critical data for the regarded twelve 2CLJDQ model fluids are 
given numerically in Table [21 which also contains the acentric factor uj 

and the critical compressibility Z,. = p*/ (p^T*). 

Fig. [9] illustrates the influence of the elongation on the enthalpy of vaporization. An 
increasing elongation reduces the enthalpy of vaporization. As can be expected (and not 
shown here), increasing polar moments increase the enthalpy of vaporization. 

Thermodynamic consistency of the simulation data was validated with the Clausius- 
Clapeyron equation 
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dlnpcr 
dT 



p„T{l/p" -l/p'Y 



Ah- 



(10) 



The vapor pressure correlation, cf. Eq. ([6]), was used for the left hand side of Eq. ( ITOl) . 
while the right hand side was calculated from present simulation data. Within statistical 
uncertainties, Eq. (|TOll is fulfilled almost throughout. 

Results from the present study can be compared with other data. Dubey and O'Shea 
[To] studied VLE of model fluids with both point dipoles and point quadrupoles embedded 
in one Lennard- Jones site, but unfortunately out of our grid of molecular parameters. In 
previous work of our group P, VLE of 2CLJQ model fluids were studied, which is a 
limiting case here. Fig. fTOl presents a comparison of the critical temperature and critical 
density for different /x*^ for a specified combination of L* = and Q*"^ = 4. It can be seen 
that present 2CLJDQ data agrees well with the other critical data (including Dubey and 
O'Shea) [HIJU]. 

5 Equation of state theory 

An EOS for 2CLJDQ fluids can be developed with the help of the VLE simulation results 
discussed above by applying a perturbation theory following the intermolecular potential, 
cf. Eqs. (11]) and ffTTl) . Using the non-polar 2CLJ fluid as a reference fluid, each of the 
polar contributions to the intermolecular potential, i.e. u^^, u^^ and u^^ , leads to a 
perturbation expansion A^^ , y^QQ ^n^j j{DQ ^.j^g Helmholtz energy A. The EOS written 
in the residual Helmholtz energy then reads 




(11) 



NksT NkeT NkeT NkeT NksT' 
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where A'^^^'^ is the residual Helmhohz energy of the 2CLJ reference fluid. The EOS of 
this reference fluid, ^^clj^ easily be constructed with Wertheim's Thermodynamic 
Perturbation Theory of first order (TPTl) [251 EZl EHl 122] • Although the TPTl initially 
resulted in a description of tangent-sphere fluids, rather than fused-sphere configurations 
(such as the 2CLJ model), the TPTl was shown to form a framework, where the 2CLJ 
fluid behavior can be recovered with high accuracy. A relation mapping the elongation L* 
of the 2CLJ fluid to the (noninteger) "number of tangent-spheres" m was earlier reported 
[9]. The other parameters of 2CLJDQ fluids, cr, e, fi* and Q*, also need to be converted 
to those of the tangent-sphere model. This conversion is straightforward and it is referred 
to Table 1 of ref. [9J. The superscript "TS" is used here to indicate parameters of the 
tangent-sphere model, e.g. it is cr = cr'^^ . The EOS of the 2CLJ fluid is then 

^^-—=m--—— + (l-m)-\ng^^(a), (12) 

where g^'^{a) is the value of the radial distribution function at the distance r = cr as 
proposed by Johnson et al. [68] . 

Expressions for A*^*^ and A^^ where reported earlier [9l |8] as third order perturbation 
expansions written in the Pade approximation. To develop an expression for A^*^, an 
analogous procedure is followed here. The EOS contribution A^^ is presented below for 
the case of mixtures, since mixtures are considered in section 6. The Helmholtz energy 
contribution A^^ to Eq. ( fTTj) is given as 

NksT l-Ag^Q/A^^' 
with A^'^ and A^'^ being the second and third order perturbation terms, respectively. 
For linear and symmetric molecules, which are considered here, the second order term 
A2^ reads in dimensionless form 
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where Xi is the mole fraction of component i, fi*^^"^ = fif / {rriiej^^ af^) denotes the dimen- 
sionless squared dipole moment in the tangent sphere framwork, and Q*^'^^ = Qf / {rriiejf a^^^ 
the respective dimensionless squared quadrupole moment. The unhke Lennard- Jones pa- 
rameters efj^ and aij were determined according to the Berthelot-Lorentz combing rules, 
with e^'^ = ieJ^e^^Y^'^ and aij = (cTj + (jj) /2. The abbreviation J^^*^ denotes inte- 
grals over angles and radius with an integrand being the angle- dependent part of the 
intermolecular potential weighted with the reference-fluid pair-correlation function. An 
analytic expression for J^^^ is not easily available and following [8], [9] , simple power func- 
tions were assumed for J^j*^, where the coefficients were adjusted to the present VLE 
simulation results. 

The third order term accounts for three-body effects and consists of two contributions 
A^^^ = A^^^ + A^^^^. Similar to J^^*^ in the second order term, both terms A^^^ and 
j^DQQ (^Qjy^g^^Yi integrals over the three-body correlation functions. The available data 
base, including the simulation results of this study, is not sufficiently broad to justify an 
independent adjustment of expressions for both of these integrals and therefore the third 
order term is simplified to 

NksT ' '''''' ksT ksT ksT af^ala^, 

f ,,*TS2,,*TS2n*TS2 , ^, , *T S2 ri*T S2 ri*T S2\ jDQ p-n 
' [H'i H'j Vfe -I- a ■ /ij i^j (Vfc ) 'J3,ijk- U^J 

It is thereby assumed that the integrals Js^ijk, which appear in both and 
have a similar density dependence and that the difference between them can be captured 
by an empirical factor a in Eq. fllSp . This is certainly a pragmatic approach and may 
need refinement once a broader data base of simulation results is available. For the 
integrals J^^*^ and J^^j^ simple power functions in density and a rudimentary temperature 
dependence in the second order term are assumed 



n=0 
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4 

•^3,!jk ~ '^^(^n,ijkV , (17) 
n=0 

where rj denotes the packing fraction, which is a dimensionless density. The relation 
between i] and p depends on the considered EOS, cf. ref. [9]. The coefficients in Eqs. ( fT6i) 
and ( fT71) depend on the chain length m with 

rriij — 1 rriij — 1 rriij — 2 

an,ij = flOn H Oi„ H a2n, (18) 

hn,ij = bon + — -bin + — -— ^&2n, (19) 

Cn,ijk = Con H Cin, (20) 

and the combining rules of the chain length are given by 



rriij 



^ {mi-rrijf^ (21) 
rriijk = {rrii ■ rrij ■ rrikf^^ (22) 

The simulation data of vapor pressure and saturated densities given in Table [U was used 
to adjust the EOS constants in Eqs. ffT^ and ffTSj) to fl^ . The contributions due to the 
quadrupole-quadrupole interactions and the dipole-dipole interactions needed for optimiz- 
ing these constants were taken from our earlier studies [HI E]- Eleven constants, namely 
a in Eq. (fT5|) as well as ao„, &on and Co„ in Eqs. ( IT8|) to (!20l) . were adjusted in a ffist 
step to the present simulation results for spherical 2CLJDQ fluids, i.e. with L* = 0. The 
remaining fifteen constants ai„, a2„ and bin, &2n and Ci„ were subsequently adjusted to 
present simulation results of elongated 2CL JDQ fluids, cf. Table [3l 

For the limiting case of non-polar 2CLJ fluids the EOS overpredicts the critical point. 
During the adjustement of the model constants, the critical point was thus not enforced. 
A comparison of the EOS to simulation data for saturated liquid and vapor denities is 
shown in Fig. [T]to|ll Apart from systematic deviations around the critical point, the EOS 
describes the data with good accuracy. A comparison of the EOS to simulation data for 
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vapor pressures is given in Fig. El to [HI Some deviations become apparant for the highest 
elongation of L* = 1, while the vapor pressure curves of spherical 2CLJDQ fluids are 
found to be well described. 



6 Application to polar mixtures 

The proposed EOS can readily be applied to mixtures of dipolar and quadrupolar flu- 
ids. To assess the capability of the proposed EOS regarding the mixed polar interac- 
tion, it was compared to simulative and experimental VLE of the following five binaries: 
C2H2+R152a, R142b+R113, RI2+CO2, R22+R142b and Propylene+R115. 

Five of these components are strongly quadrupolar C2H2 {Q = 5.1 DA), R113 (13.0 
DA), Propylene (5.9 DA), R115 (9.2 DA) and CO2 (3.8 DA) and were modelled by 2CLJQ 
models in [21 H] . The numbers in parentheses indicate the polar moment of the molecular 
model, respectively. The remaining four components R152a (/i = 2.7 D), R142b (3.0 
D), R22 (2.3 D) and R12 (2.3 D) are dipolar and were modelled by 2CLJD models in 
[3]. Hence, the first three mixtures mentioned above exhibit the dipole-quadrupole cross- 
interaction, whereas the remaining two are dipolar or quadrupolar only and are limiting 
cases here. Simulative VLE calculations were performed on the basis of these molecular 
models in [69], where one state independent binary parameter ^ for the unlike Lennard- 
Jones interaction 



was adjusted to one experimental vapor pressure following the procedure proposed in 0. 
Table 4 contains the binary parameters derived in [02] and Table 5 compiles the binary 
VLE data from simulation in numerical form. As can be seen in Figs. 11 to 15, the 
simulative approach agrees excellently to the experimental data which exhibits in part a 
qualitatively different phase behavior. 




(23) 
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The adopted molecular model (2CLJDQ) is of course simple for the considered real 
fluids; the most delicate assumptions are: Firstly, the polarizability is not accounted for 
and vacuum-value dipole moments are used. Secondly, the multipole moments are not 
represented by (distributed) charges but by point-multipoles. Thirdly, the multipolar 
moments are assumed to be aligned along the molecular axis. The good agreement of the 
simulation results to the experimental phase behavior however suggests that important 
characteristics of the real fluids are captured. 

To validate the EOS containing the new dipole-quadrupole contribution, these molec- 
ular models were used with exactly the same parameters without any further adjustment. 
The results are also shown in Figs. 11 to 15. It can be seen that the EOS is in excellent 
agreement with the simulations, usually within the statistical uncertainty of the simula- 
tion data which is often within symbol size. Fig. 12 also gives a comparison with the 
EOS suggested by Weingerl and Fischer [25]- This model (dotted line) is found in good 
but not entirely quantitative agreement. 

The components of, e.g., the mixture R142b+R113 have high dipole and quadrupole 
moments, leading to considerable dipole-quadrupole cross-interactions in the mixture. 
The observation that the EOS predicts the mixture phase behavior well suggests that the 
polar cross-interactions are adequately captured by the new EOS model. This conclusion 
is supported by a calculation also shown in Fig. 12, where the dipole-quadrupole cross- 
iteractions between the two compounds are set to zero. This is done in order to illustrate 
the contribution that is due to the here proposed expression for A^*^. The EOS is for 
that case (dashed line) not in the vicinity of the data. 
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7 Conclusion 



The VLE of dipolar and quadrupolar two-center Lennard- Jones fluids was studied by 
molecular simulation to derive a new EOS contribution for the dipole-quadrupole inter- 
actions of non-spherical fluids. The vapor pressure, saturated densities and saturated 
enthalpies were determined with the NpT -\- test particle method applying a gradual in- 
sertion scheme in cases of high polar densities. The EOS expression has the form of a 
perturbation theory of third order where coefficients were adjusted to the bulk properties 
of the fluids considered in this work. The EOS was subsequently applied to mixtures 
without adjustable parameters where an excellent agreement to experiment and simula- 
tion data was found. 
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List of Symbols 



Latin alphabet 



a 


equation of state parameter 


a 


interaction site index 


A 


Helmholtz energy 


b 


equation of state parameter 


b 


interaction site index 


c 


equation of state parameter 


c 


coefficient of correlation function 


C 


coefficient of correlation function 


h 


enthalpy 


i 


molecule index 


3 


molecule index 


J 


integral over intermolecular potential 


ks 


Boltzmann constant 


L 


molecular elongation 


m 


chain length 


N 


number of particles 


P 


pressure 


Q 


quadrupolar moment 


r 


site-site distance 


rc 


center-center cut-off radius 


T 


temperature 


u 


pair potential 


X 


mole fraction in the saturated liquid 


y 


mole fraction in the saturated vapor 


z 


compressibility 


Vector properties 


r 


position vector 




orientation vector 
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Greek alphabet 



a equation of state parameter 

Ahv enthalpy of vaporization 

At integration time step 

e Lennard- Jones energy parameter 

eg relative permittivity of dielectric continuum 

T) packing fraction 

9 angle of nutation 

II dipolar moment 

p density 

a Lennard- Jones size parameter 

(pij azimuthal angle between the axis of molecules i and j 

Lo acentric factor 

Subscript 

c property at critical point 

D dipole 

i component index 

J component index 

k component index 

Q quadrupole 

o vapor-liquid coexistence 

2CLJ two-center Lcnnard-Jones 

2CLJD two-center Lennard- Jones plus point dipole 

2CLJQ two-center Lennard- Jones plus point quadrupole 

2CLJDQ two-center Lennard- Jones plus point dipole and point quadrupole 



Superscript 



* reduced property 

' saturated liquid 

" saturated vapor 

id ideal gas property 

res residual property 

TS tangent sphere 
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Table 1: Vapor-liquid equilibrium data of twelve 2CLJDQ model fluids. The number 
parentheses indicates the statistical uncertainty in the last decimal digit. 











p"* 










L* = 0, = 


6, g 


*2 = 2 














3.78268 


0.0059 


(6) 


0.8943 (2) 


0.0013 


(3) 


-42.01 


(1) 


0.2 


(2) 


4.12656 


0.0152 


(9) 


0.8510 (2) 


0.0038 


(3) 


-40.38 


(1) 


-0.5 


(2) 


4.47044 


0.032 


(1) 


0.8240 (3) 


0.0079 


(4) 


-38.74 


(1) 


-1.3 


(2) 


4.81432 


0.059 


(2) 


0.7870 (3) 


0.0138 


(6) 


-37.10 


(1) 


-1.8 


(1) 


5.15820 


0.104 


(2) 


0.7463 (4) 


0.0237 


(6) 


-35.35 


(2) 


-2.60 


(7) 


5.50208 


0.170 


(2) 


0.6997 (4) 


0.0386 


(5) 


-33.42 


(2) 


-3.85 


(7) 


5.84596 


0.260 


(2) 


0.6474 (8) 


0.0605 


(7) 


-31.29 


(3) 


-5.41 


(9) 


6.18984 


0.381 


(3) 


0.584 (1) 


0.097 


(2) 


-28.74 


(5) 


-7.8 


(2) 


6.31849 


0.437 


(4) 


0.552 (2) 


0.116 


(3) 


-27.51 


(7) 


-8.9 


(2) 


L* = 0, 


fi*'' = 


6, g 


*2 = 4 














4.18204 


0.005 


(2) 


0.96057 (3) 


0.0013 


(4) 


-54.11 


(2) 


-0.5 


(1) 


4.59192 


0.008 


(2) 


0.92333 (3) 


0.0018 


(5) 


-51.84 


(2) 


-0.5 


(1) 


4.97458 


0.023 


(2) 


0.88729 (2) 


0.0050 


(3) 


-49.70 


(1) 


-1.18 (3) 


5.35724 


0.044 


(3) 


0.84937 (3) 


0.0090 


(7) 


-47.57 


(1) 


-1.8 


(2) 


5.73990 


0.080 


(3) 


0.80792 (4) 


0.0157 


(7) 


-45.32 


(2) 


-2.6 


(1) 


6.12256 


0.137 


(3) 


0.76301 (4) 


0.0277 


(9) 


-42.97 


(2) 


-4.2 


(1) 


6.50522 


0.227 


(3) 


0.71269 (6) 


0.0469 


(1) 


-40.42 


(3) 


-5.8 


(2) 


6.84334 


0.320 


(6) 


0.6582 (2) 


0.067 


(2) 


-37.75 


(8) 


-7.6 


(2) 


7.03343 


0.394 


(5) 


0.6257 (2) 


0.085 


(3) 


-36.25 


(5) 


-8.9 


(3) 


7.22352 


0.477 


(5) 


0.5852 (2) 


0.112 


(3) 


-34.33 


(9) 


-11.0 


(3) 


L* = 0, = 


12, g*2 = 2 














4.36466 


0.0049 


(5) 


0.9208 (3) 


0.0012 


(1) 


-53.68 


(1) 


-0.62 (6) 


4.76145 


0.007 


(1) 


0.8868 (3) 


0.0015 


(2) 


-51.82 


(1) 


-0.6 


(1) 


5.32623 


0.029 


(2) 


0.8367 (3) 


0.0060 


(4) 


-49.18 


(1) 


-1.8 


(2) 


5.73594 


0.052 


(2) 


0.7971 (3) 


0.0102 


(4) 


-47.19 


(2) 


-2.51 (8) 


6.14565 


0.096 


(2) 


0.7556 (4) 


0.0188 


(5) 


-45.15 


(2) 


-3.9 


(1) 


6.55536 


0.161 


(3) 


0.7088 (5) 


0.0318 


(6) 


-42.91 


(2) 


-5.6 


(1) 


6.96507 


0.249 


(3) 


0.6550 (7) 


0.0509 


(9) 


-40.41 


(3) 


-7.7 


(1) 


7.37478 


0.380 


(4) 


0.590 (1) 


0.087 


(2) 


-37.42 


(4) 


-11.2 


(2) 


7.53896 


0.432 


(4) 


0.554 (2) 


0.100 


(2) 


-35.81 


(8) 


-12.1 


(2) 


7.57963 


0.445 


(5) 


0.544 (2) 


0.106 


(3) 


-35.40 


(7) 


-12.7 


(3) 
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Table [2 continued. 





Pi 


P 






p * 




^/rcs 




^//rcs* 


L* = 0, = 


Q*"^ = 4 
















5.25366 


0.003 (3) 


0.9433 


(3) 





0006 


(5) 


-63.53 


(2) 


-0.4 (3) 


5.69146 


0.015 (2) 


0.9075 


(3) 





0027 


(4) 


-61.18 


(2) 


-1.4 (2) 


6.12927 


0.030 (4) 


0.8696 


(3) 





0053 


(9) 


-58.80 


(1) 


-2.1 (5) 


6.56707 


0.057 (3) 


0.8296 


(3) 





0099 


(7) 


-56.38 


(2) 


-3.3 (2) 


7.00488 


0.103 (6) 


0.7858 


(5) 





018 


(1) 


-53.81 


(2) 


-4.5 (4) 


7.44268 


0.185 (4) 


0.7379 


(6) 





034 


(1) 


-51.08 


(3) 


-7.5 (3) 


7.88049 


0.290 (5) 


0.6834 


(9) 





049 


(2) 


-48.07 


(4) 


-10.1 (2) 


8.09939 


0.347 (6) 


0.651 


(1) 





060 


(4) 


-46.33 


(4) 


-10.9 (6) 


8.45000 


0.478 (8) 


0.591 


(2) 





104 


(4) 


-43.21 


(7) 


-15.2 (5) 


L* = 0.505, fi*^ = 


6, Q*2 -_ 


= 2 














1.94604 


0.0016 (3) 


0.5765 


(3) 





0008 


(2) 


-23.96 


(1) 


-0.21 (5) 


2.12295 


0.0028 (5) 


0.5542 


(3) 





0014 


(2) 


-22.96 


(1) 


-0.30 (4) 


2.29986 


0.0064 (7) 


0.5313 


(2) 





0030 


(3) 


-21.979 


(9) 


-0.48 (6) 


2.47678 


0.0142 (8) 


0.5076 


(2) 





0062 


(4) 


-21.011 


(9) 


-0.80 (5) 


2.65369 


0.026 (2) 


0.4816 


(6) 





0113 


(9) 


-19.99 


(2) 


-1.3 (1) 


2.83060 


0.047 (2) 


0.4525 


(4) 





021 


(1) 


-18.88 


(1) 


-2.2 (3) 


3.00752 


0.074 (1) 


0.4184 


(6) 





032 


(1) 


-17.64 


(2) 


-2.7 (2) 


3.18443 


0.110 (1) 


0.377 


(1) 





050 


(1) 


-16.20 


(4) 


-4.0 (2) 


3.27289 


0.131 (1) 


0.351 


(4) 





064 


(2) 


-15.29 


(9) 


-4.7 (1) 


L* = 0.505, fi*'' = 


6, Q*'' -- 


= 4 














2.30000 


0.0011 (5) 


0.5893 


(2) 





0005 


(2) 


-28.78 


(1) 


-0.2 (1) 


2.49167 


0.0056 (7) 


0.5664 


(2) 





0024 


(3) 


-27.53 


(1) 


-0.67 (7) 


2.68333 


0.009 (2) 


0.5426 


(3) 





0037 


(9) 


-26.30 


(1) 


-0.83 (2) 


2.87500 


0.017 (3) 


0.5162 


(3) 





007 


(1) 


-25.01 


(1) 


-1.2 (2) 


3.06667 


0.033 (4) 


0.4884 


(4) 





013 


(2) 


-23.69 


(1) 


-1.9 (3) 


3.25833 


0.058 (2) 


0.4562 


(4) 





023 


(1) 


-22.25 


(2) 


-2.8 (1) 


3.45000 


0.093 (2) 


0.4198 


(7) 





039 


(2) 


-20.66 


(2) 


-4.3 (2) 


3.54583 


0.117 (3) 


0.3983 


(9) 





052 


(2) 


-19.79 


(3) 


-5.1 (2) 


3.64167 


0.143 (2) 


0.372 


(1) 





070 


(2) 


-18.77 


(4) 


-6.4 (1) 
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Table [U continued. 





P*a 


P'* 




p"* 






j^/r res" 


L* = 0.505, fi*'^ = 


12, Q*'' -- 


= 2 












2.23188 


0.0011 (3) 


0.5917 


(2) 


0.0005 (2) 


-30.816 


(9) 


-0.63 


(9) 


2.44569 


0.0016 (6) 


0.5690 


(2) 


0.0007 (3) 


-29.62 


(1) 


-0.3 


(3) 


2.64949 


0.0040 (6) 


0.5466 


(2) 


0.0016 (3) 


-28.47 


(1) 


-0.74 


(9) 


2.85330 


0.0093 (9) 


0.5226 


(3) 


0.0036 (4) 


-27.33 


(1) 


-1.2 


(2) 


3.06000 


0.0189 (7) 


0.4964 


(3) 


0.0072 (3) 


-26.10 


(1) 


-1.88 


(9) 


3.26092 


0.0304 (2) 


0.4687 


(4) 


0.012 (1) 


-24.87 


(2) 


-2.9 


(3) 

V'-'/ 


3.46472 


0.055 (1) 


0.4384 


(5) 


0.0213 (6) 


-23.57 


(2) 


-3.76 


(9) 


3.66853 


0.080 (2) 


0.399 


(1) 


0.034 (1) 


-21.97 


(3) 


-5.0 


(1) 


3.85506 


0.116 (3) 


0.354 


(2) 


0.050 (3) 


-20.20 


(5) 


-6.4 


(3) 


L* = 0.505, jj*'^ = 


12, g*^ : 


= 4 












2.47462 


0.00052 (1) 


0.6181 


(2) 


0.0001 (2) 


-37.11 


(1) 


-0.5 


(4) 


2.69958 


0.00124 (1) 


0.5958 


(2) 


0.00034 (5) 


-35.67 


(1) 


-0.1 


(2) 


2.92455 


0.00322 (1) 


0.5722 


(3) 


0.00110 (3) 


-34.23 


(1) 


-1.17 


(6) 


3.14951 


0.00787 (2) 


0.5477 


(3) 


0.00275 (1) 


-32.79 


(1) 


-1.54 


(2) 


3.37448 


0.01649 (8) 


0.5214 


(3) 


0.00562 (5) 


-31.32 


(1) 


-2.22 


(3) 


3.59944 


0.0304 (1) 


0.4926 


(4) 


0.01021 (7) 


-29.82 


(2) 


-3.11 


(4) 


3.82441 


0.049 (3) 


0.4593 


(6) 


0.0177 (6) 


-28.10 


(3) 


-4.5 


(1) 


4.04937 


0.080 (4) 


0.422 


(1) 


0.029 (3) 


-26.29 


(4) 


-5.8 


(5) 


4.16186 


0.105 (4) 


0.400 


(1) 


0.042 (3) 


-25.28 


(4) 


-7.2 


(5) 


4.27434 


0.132 (3) 


0.371 


(2) 


0.056 (2) 


-23.99 


(6) 


-8.4 


(3) 


L* = 1, 


= 6, g 


*2 = 2 














1.35000 


0.00013 (1) 


0.4694 


(2) 


0.00009 (3) 


-23.03 


(2) 


-0.25 


(1) 


1.48000 


0.00060 (1) 


0.4466 


(3) 


0.00035 (4) 


-21.50 


(2) 


-0.29 


(1) 


1.60000 


0.00187 (3) 


0.4257 


(3) 


0.00126 (5) 


-20.22 


(2) 


-0.56 


(6) 


1.72000 


0.00448 (2) 


0.4036 


(3) 


0.00280 (6) 


-18.98 


(1) 


-0.93 


(4) 


1.85000 


0.00989 (3) 


0.3781 


(2) 


0.00629 (3) 


-17.63 


(1) 


-1.40 


(2) 


1.96000 


0.01775 (7) 


0.3549 


(3) 


0.01138 (6) 


-16.50 


(1) 


-1.93 


(1) 


2.09000 


0.0326 (2) 


0.3233 


(5) 


0.0218 (2) 


-15.06 


(2) 


-2.80 


(4) 


2.21000 


0.0522 (3) 


0.2840 


(9) 


0.0370 (5) 


-13.46 


(3) 


-3.86 


(8) 


2.26000 


0.0627 (4) 


0.251 


(2) 


0.0456 (7) 


-12.35 


(5) 


-4.34 


(9) 
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Table [TJ continued. 





Pi 


P'* 




p"* 










L* = 1, 


p*^ = 6, Q*' 


= 4 
















1.84000 


0.0005 (1) 


0.4787 


(8) 


0.00033 


(2) 


-30.34 


(8) 


-1.0 


(2) 


2.00000 


0.0019 (1) 


0.4505 


(3) 


0.00100 


(2) 


-28.14 


(2) 


-1.38 (9) 


2.13000 


0.0047 (4) 


0.4280 


(4) 


0.00263 


(4) 


-26.43 


(3) 


-2.12 (9) 


2.25000 


0.0101 (6) 


0.4052 


(4) 


0.0055 


(4) 


-24.85 


(3) 


-2.9 


(2) 


2.40000 


0.0238 (3) 


0.3731 


(4) 


0.0136 


(4) 


-22.79 


(2) 


-4.4 


(2) 


2.55000 


0.040 (1) 


0.3367 


(6) 


0.0236 


(1) 


-20.68 


(3) 


-5.1 


(2) 


2.60000 


0.049 (1) 


0.3251 


(8) 


0.0305 


(3) 


-19.94 


(3) 


-5.8 


(3) 


2.68000 


0.065 (5) 


0.300 


(1) 


0.047 


(8) 


-18.70 


(4) 


-7.1 


(8) 


2.75000 


0.091 (8) 


0.267 


(2) 


0.07 


(1) 


-17.13 


(7) 


-9 


(1) 


L* = 1, 


/i*' = 12, Q 


*2 = 2 
















1.64000 


0.0001 (1) 


0.4773 


(5) 


0.00006 


(1) 


-30.49 


(4) 


-1.1 


(3) 


1.80000 


0.0004 (1) 


0.4572 


(3) 


0.00018 


(2) 


-29.02 


(2) 


-1.35 


(9) 


1.94000 


0.0011 (1) 


0.4350 


(3) 


0.00065 


(1) 


-27.31 


(3) 


-1.76 


(5) 


2.10000 


0.0037 (3) 


0.4096 


(3) 


0.0014 


(3) 


-25.52 


(2) 


-2.0 


(4) 


2.25000 


0.0088 (5) 


0.3820 


(3) 


0.0046 


(6) 


-23.75 


(2) 


-3.2 


(3) 


2.37000 


0.0148 (1) 


0.3588 


(4) 


0.00853 


(9) 


-22.36 


(2) 


-4.03 


(4) 


2.50000 


0.0276 (5) 


0.3305 


(5) 


0.0173 


(5) 


-20.77 


(2) 


-5.3 


(1) 


2.70000 


0.0543 (5) 


0.264 


(2) 


0.0353 


(7) 


-17.64 


(6) 


-6.9 


(1) 


2.80000 


0.073 (1) 


0.228 


(3) 


0.054 


(1) 


-16.01 


(9) 


-7.9 


(1) 


L* = 1, 


/i*' = 12, Q 


*2 = 4 
















2.14000 


0.00010 (1) 


0.4897 


(9) 


0.00008 


(1) 


-40.21 


(7) 


-3.2 


(3) 


2.31000 


0.00051 (1) 


0.4790 


(7) 


0.00032 


(1) 


-39.43 


(7) 


-4.2 


(2) 


2.48000 


0.00157 (2) 


0.461 


(1) 


0.00082 


(3) 


-38.0 


(1) 


-4.4 


(3) 


2.61000 


0.00393 (7) 


0.4338 


(6) 


0.00191 


(3) 


-35.64 


(4) 


-5.56 (5) 


2.75000 


0.0081 (2) 


0.4090 


(6) 


0.0041 


(1) 


-33.55 


(4) 


-6.8 


(1) 


2.97000 


0.023 (1) 


0.3646 


(6) 


0.0125 


(9) 


-30.24 


(3) 


-8.8 


(3) 


3.08000 


0.032 (1) 


0.3411 


(8) 


0.0177 


(9) 


-28.60 


(4) 


-9.4 


(2) 


3.19000 


0.047 (2) 


0.3132 


(9) 


0.027 


(2) 


-26.79 


(4) 


-10.3 


(4) 
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Table 2: Critical temperature, density, pressure and compressibility as well as acentric 
factor of twelve 2CLJDQ model fluids. 





L* 









0.505 


1 












6 


12 


6 


12 


6 


12 








6.651 


7.937 


3.445 


4.067 


2.344 


2.863 


C 






0.319 


0.3128 


0.1982 


0.1929 


0.1508 


0.1350 


Pc 




2 


0.5995 


0.6125 


0.1844 


0.1692 


0.0822 


0.0897 


Pc 






0.2823 


0.2467 


0.2701 


0.2157 


0.2324 


0.2320 


Zc 






0.1254 


0.1733 


0.1987 


0.2423 


2.4890 


2.8628 


OJ 






7.604 


8.944 


3.832 


4.501 


2.861 


3.424 


C 






0.3288 


0.3229 


0.2086 


0.1989 


0.1532 


0.1481 


Pi 




4 


0.6797 


0.7036 


0.2121 


0.2382 


0.1335 


0.1112 


Pi 






0.2718 


0.2436 


0.2654 


0.2660 


0.3045 


0.2191 








0.2057 


0.2783 


0.3296 


2.4341 


3.2371 


3.8008 
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Table 3: Model constants of the new dipole-quadrupole EOS contribution. In Eq. (|TH|) . 
the parameter a=l. 19374. 



i 




au 


a2i 


boi 


bit 


b2i 


Coi 


Cli 





0.6970950 


-0.6734593 


0.6703408 


-0.4840383 


0.6765101 


-1.1675601 


7.846431 


-20.72202 


1 


-0.6335541 


-1.4258991 


-4.3384718 


1.9704055 


-3.0138675 


2.1348843 


33.42700 


-58.63904 


2 


2.9455090 


4.1944139 


7.2341684 


-2.1185727 


0.4674266 





4.689111 


-1.764887 


3 


-1.4670273 


1.0266216 
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Table 4: Binary interaction parameters taken from [U^ . 



Mixture 




C2H2+R152a 


1.090 


R142b+R113 


0.952 


CO2+RI2 


0.927 


R22+R142b 


0.985 


Propylene+R115 


0.948 
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Table 5: Vapor-liquid equilibrium data for five binary mixtures taken from Pure 
substance vapor pressures were obtained via correlations provided in The number 

in parentheses indicates the statistical uncertainty in the last decimal digit. 



Xi / mol/mol 


p 1 MPa 


y\ 1 mol/mol 


C2H2+R152a ( 


1+2) at 303.2 K 





0.69 





0.128 


0.87 (5) 


0.32 (2) 


0.424 


1.77 (5) 


0.76 (2) 


0.569 


2.45 (8) 


0.88 (2) 


0.801 


4.01 (6) 


0.950 (4) 


1 


5.62 


1 


R142b+R113 ( 


1+2) at 373 K 





0.45 





0.252 


0.90 (3) 


0.57 (2) 


0.502 


1.27 (4) 


0.77 (2) 


0.751 


1.63 (4) 


0.890 (5) 


1 


2.00 


1 


CO2+RI2 (1+2) at 273 K 





0.30 





0.147 


0.86 (2) 


0.65 (1) 


0.388 


1.70 (2) 


0.849 (6) 


0.550 


2.20 (2) 


0.899 (4) 


0.714 


2.67 (2) 


0.932 (4) 


1 


3.52 


1 


R22+R142b (1+2) at 328 


15 K 





0.75 





0.218 


1.02 (1) 


0.40 (1) 


0.455 


1.36 (1) 


0.66 (1) 


0.560 


1.50 (3) 


0.730 (8) 


0.661 


1.66 (2) 


0.806 (5) 


0.806 


1.86 (2) 


0.900 (3) 


1 


2.17 


1 


Propylene+R115 (1+2) at 298 K 





0.95 





0.113 


0.98 (4) 


0.175 (7) 


0.316 


1.18 (3) 


0.386 (8) 


0.549 


1.24 (2) 


0.59 (1) 


0.788 


1.26 (2) 


0.780 (7) 


1 


1.19 


1 
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Figure 1: Vrabec and Gross 
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Figure 3: Vrabec and Gross 
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Figure 4: Vrabec and Gross 
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Figure 5: Vrabec and Gross 
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Figure 7: Vrabec and Gross 
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Figure 9: Vrabec and Gross 
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Figure 10: Vrabec and Gross 
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Figure 12: Vrabec and Gross 
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Figure 13: Vrabec and Gross 
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Figure 14: Vrabec and Gross 
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Figure 15: Vrabec and Gross 
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